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Abstract 

We generalize some widely used mother wavelets by means of the 
g-exponential function e? = [1 + (1 — q)x] 1 ^ 1 ~ q ^ (q 6 M, ef = e x ) that 
emerges from nonextensive statistical mechanics. Particularly, we de- 
fine extended versions of the mexican hat and the Morlet wavelets. We 
also introduce new wavelets that are ^-generalizations of the trigono- 
metric functions. All cases reduce to the usual ones as q — > 1. Within 
nonextensive statistical mechanics, departures from unity of the en- 
tropic index q are expected in the presence of long-range interactions, 
long-term memory, multi- fractal structures, among others. Consis- 
tently the analysis of signals associated with such features is hopefully 
improved by proper tuning of the value of q. We exemplify with the 
WTMM Method for mono- and multi- fractal self-affine signals. 

PACS: 02.30.-f, 02.30.Nw, 05.20.-y, 05.45.Tp, 
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1 Introduction 



The basic idea behind the analysis of a temporal or spatial signal by a Fourier 
transform is similarity. For instance, the inner product (<f) n (x), f(x)) ex- 
presses how similar are the function f(x) and the n th element of the basis 
{4> n {x)} (the kernel of the transform). The kernel of Fourier transform is a 
plane wave and here lies its simplicity, but also its limitation. Rigorously 
speaking, a plane wave exists everywhere in the universe and at all (past and 
future) times. Hence there are no plane waves in nature. All physical sig- 
nals are temporary and spatially limited. As a consequence of the infinitely 
extended kernel, the Fourier transform is unable to satisfactorily determine 
when a burst has occured, or where edges of images are located, for example 
(the Gibbs phenomenon, see, for instance, PUI2])- 

The first attempt to overcome limitations of Fourier analysis was that 
of Gabor j3J, who introduced a windowed Fourier transform |4.. He used a 
modulation function g(u — t) in order to localize the signal in a time interval 
[t — T,t] (supp g C [— T, 0]). Translations in time would cover the whole 
signal. With this, he achieved time-frequency localization. But a problem 
still remains with the use of a fixed scale window: signal details much smaller 
than the window width T are detected, but not localized. They appear in 
the frequency behaviour of the windowed Fourier transform in a similar way 
it would appear in the usual Fourier transform. Signal features much larger 
than T, on the other hand, appear in the time behaviour of the windowed 
Fourier transform (i.e., they are not detected). 

To avoid this problem, it is sufficient a scale free transform. Wavelet anal- 
ysis was developed to accomplish this goal. The windowed Fourier transform 
uses a fixed size window and fills it with oscillations of different frequencies 
(consequently, varies the number of oscillations within the window). The 
wavelet transform uses a function with fixed number of oscillations and vary 
the width of the window. Dilations and translations generate a complete 
analysis of the signal. This procedure automatically uses small windows to 
identify high-frequency components of a signal, and large windows for low- 
frequency components. Wavelet analysis provides a time-scale (instead of a 
time-frequency) localization. 

Wavelet analysis was developed more than a century and a half after the 
pioneering work of Fourier j^J. Since the 1980's it has grown so fast that it 
has already consolidated as a field of research of its own. 
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Another emerging science is nonextensive statistical mechanics, that also 
began in the 1980's, with the formulation of the concept of a nonextensive 
entropy jH], that generalizes the Boltzmann-Gibbs entropy, as a basis of a 
generalization of standard statistical mechanics itself. Here, it took also 
more than a century since the first formulation of the concept of entropy. 
Nonextensive statistical mechanics (HI El M E] seems to give an underlying 
formal basis for a variety of complex phenomena, such as anomalous diffusion 
PHI ED E2 E31 EH E3 EE], self-gravitating systems H3 EH, turbulence [HI 
EDI Ell 122 E31 ES ES] cosmic rays [23 EH, among others (for un updated 
bibliography, see [28]). 

The definition of the nonextensive g-entropy jU] is 

-i v-^iy q 

S = k ZA^h (fc > o). (l) 
q — 1 

The entropic index g characterizes the generalization and the Boltzmann- 
Gibbs-Shanon entropy, Si = —kJ^iPi^Pi, is recovered at q — > 1. The 
canonical ensemble is obtained by maximizing with the norm constraint 
(5^i=i P« = 1) an d also imposing the constancy of the normalized g-expectation 
value of the energy [2:9 

w 

U q = J2^"\ (2) 

where {ej}^ are the eigenvalues of the Hamiltonian of the system, and 
x are the escort probabilities, defined by [SHj 



pW) — ft (o\ 



,(«) = _P_ 

2^=i Pj 

with pi being the probability associated with the i th state, given by 



Pi = (4) 

If q < 1, pi = whenever [1 — (1 — q)/3'ei\ < (cut-off condition). The 
partition function Z' q is defined as 

w 

% = (1-9)/?*]^ (5) 



3 



and 

22j=iPj + 0--Q)P U 9 

where (3 is the Lagrange parameter associated with the constraint (J2J). The 
striking features of (J3J) are its power-law for q > 1, and the cut-off for 
g < 1, which generates quite different behaviours when compared with the 
Boltzmann-Gibbs exponential distribution. The meaning of the word nonex- 
tensive can be easily understood if we consider a system composed of two in- 
dependent subsystems (in the sense of probability theory, Pif +B ^ = p\ p\ i )• 
According to (JT]), the g-entropy of the composite system is given by 

g(A + B) = S (A) + S (B) + (1 _ q)S (A) S (B)_ (7) 

Clearly g / 1 yields nonadditivity, or nonextensivity. A system with q > 1 
is said to be subadditive, and if q < 1, it is superadditive. 

The functional forms of (0) and inspired the definition of generaliza- 
tions of the logarithm and exponential functions (SUE2] : 

xl ~ q ~ 1 

ln i x = 1 _ > ( 8 ) 




-q)x]i-«, if [1 + (1 - q)x) > (g) 
otherwise. 

It follows immediately that ln g (eg) = e^ 9 ^ = x. The original logarithm 
and exponential functions are recovered as the particular cases lnx x and 
ef. Many formal developments have been made concerning such functions. 
Among them we point out the generalization [33] of the celebrated Shannon's 
theorem. Obviously there are infinitely many ways of generalizing a given 
function. One commonly form of generalizing the exponential is exp ? (x) = 

En^/W? 1 - With ( n )q- = YYj=lU)q and (j)<? = ~ 1) I ' (l ~ X ) and als ° 

(0) g ! = 1. This particular expression have applications in quantum groups 
|34j . In this paper we won't explore this possibility, but rather that one given 
by Equation © above, that naturally emerges from nonextensive statistical 
mechanics. Relations between wavelet analysis and nonextensive statistical 
mechanics have already been reported, including applications in biophysics, 
as the analysis of EEG signals [23 EE1 E3 EH1 EH • 
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The purpose of this paper is to extend the use of such g-deformed func- 
tions into wavelet analysis by generalizing some widely used wavelet func- 
tions. Particularly, we generalize the mexican hat (Section |2J) and the modu- 
lated Gaussian (Section|3J). Moreover, we introduce in SectionHJa pair of even 
and odd wavelets based on the generalization of trigonometric functions. Fi- 
nally (Section|5J) we exemplify, with the Wavelet Transform Modulus-Maxima 
Method for mono- and multi-fractal self-affine signals, the possible use of the 
introduced functions, and, lastly, we present our final remarks (Section EJ. 



2 (/-Mexican Hat 

Before we discuss the properties of the wavelets based on definitions (jHJ) and 
(Q, let us indicate some basic facts about continuous ^-wavelet transform of 
a function /, T g (a, b)f, which is defined by 

/°° t h 
dxf(x)g( ). (10) 
-oo a 

In order that a given function g G L 2 (IR) can be considered as a mother 
wavelet, it must satisfy the so-called admissibility condition, that can be 
expressed in terms of its Fourier transform g{oS) by 

/oo 
du\g(u)\ 2 /\u\ < oo. (11) 
-oo 

It essentially means that g(x) has zero mean, i.e., 

/oo 
dxg(x) = (12) 
-oo 

and ensures that the the original function f(x) can be recovered from its 
wavelet transform as 

/OO /■OO 1 r — h 
db da-T g (a,b)fg( -). (13) 
-oo JO a a 

This expression shows that, if a mother wavelet g satisfies (jllj) . it generates, 
by its translations and dilations, a basis in the Hilbert space, so that any 
function / G L 2 (R) can be expressed in terms of its g— wavelet components. 
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It is well known that orthonormal wavelet bases can only be denned within 
the discrete formulation, i.e., when b G Z and a = a J , with a > and j G Z. 
So, orthonormal bases can not be defined with the help of the functions we 
discuss below. 

A simple and quite common example of continuous wavelet is the mexican 
hat (see, for instance, jJUJ EU E2|), that is generated from the Gaussian 
distribution: 

d 2 e- x ' 2 / 2 

= A(l-x 2 )e- x2/2 . (14) 

The normalization constant is given by A = 2/(tt 1 / 4 v / 3). Generalizat ion 
of Gaussian distribution within nonextensive scenario we are focusing here 
has already been made [T3J [TH]. The q- Gaussian (oc e~^ x ) unifies a great 
variety of different distributions into a single family, parameterized by q (see 
[TH| ITT)] for details): Gaussian distribution is recovered, of course, for q — 1. 
Moreover, for q = 2, we have the Cauchy-Lorentz distribution. Besides, 
q — > 3 yields a completely flat distribution, and q — > — oo yields Dirac's 5. 
This unifying character of the g-Gaussian, as well as its empirically observed 
occurence in a variety of complex phenomena (as cited in the Introdutcion) 
stimulate us to explore its potential use in wavelet analysis. One of the more 
simple applications that immediately come to our minds is to use it in order 
to generalize the mexican hat, a function that is proportional to the second 
derivative of a Gaussian. We could simply replace an ordinary Gaussian by 
a g-Gaussian, which is a legitimate choice. Instead of doing that, we prefer 
to adopt a variation of it, and take the second derivative of a power of a 
g-Gaussian, according to the recipe 

^2 |" e -y9a; a j2-g 

M x ) oc ^ , (15) 

which yields the expression for the g-mexican hat 

%{x) = A q [1 - (3 - q)(3x 2 } [e-^ 2 ] " . (16) 

With this choice we arrive, according to the expression above, at a g-Gaussian 
raised to the power g, resembling the escort probabilities, Equation (JHJ), 
widely used in nonextensive formalism [201 • We call the attention of the 
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reader that [e^] b ^ [e^ 6 ], except, of course, in the particular case q = 1. Be- 
cause of this, the semigroup property is not satisfied, but it is exactly this 
that generates its nonextensive behaviour Equation (JJJ). For some properties 
of the g-exponential and the q- logarithm, see [I3J HH EE EE HZ] • 

If — 1 < q < 1 and \x\ = [(1 — q)P]^ 1 ^ 2 , the cut-off of the g-exponential (see 



Equation Q) imposes ip q 



0. It is easily proven that, for — 1 < q < 3, 



Equation (jllj) holds, so that any ip q {x) is a well defined mother wavelet. For 
q < —1 or q > 3, ip q {x) does not satisfy the admissibility condition and can 
not be used to define a useful mother wavelet. 
The normalization constant is given by 



(3^_ 



(?-i) 5/2 r(^ 



n 1/2 



(17) 



if 1 < q < 3, and 



01/4 ( 5 _ ? )l/2 (3 + g )l/2 



7rV4^3 



:i-^) 1/2 rfe + 



r(^ + i 



1/2 



(18) 



if — 1 < q < 1. 

The function ip q (x) satisfies the admissibility condition and consistently 
recovers the mexican hat, lim^! ip q {x) = ipi{x). The range of admissible 
values for q (—1 < q < 3) is divided in three regions. For 1 < q < 3, ip q (x) is 
infinitely supported and presents a power-law tail ~ —l/\x\ 2 ^ q ^ 1 \ in marked 
contrast with the exponential tail of the original mexican hat. When q < 1, 
a cut-off naturally appears at \x c \ = [(1 — q)PY 1 ^ 2 - In the range < q < 1, 
i/j q (x c ) = 0, and when — 1 < q < 0, i/j q (x c ) diverges. For q — > —1, ^ q (x) 
coincides with the abscissa axis, except at the cut-off positions, where it 
diverges. These features introduce significant differences from the original 
mexican hat wavelet. Figure illustrates ip q {x) with (3 = 1/2. 

The Fourier transform, 



F[f{x)-y\=F{y) 



e ixy f(x) dx, 



(19) 



of i( q (x) may be found by considering (|15p. and taking into account the 



property of the Fourier transform of derivatives, T [f^ n '; y\ 



-iy) n F(y), 
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together with the Fourier transform of a g-Gaussian (see Equations 3.384 9 
and 3.387 2 of @E]). We find for 1 < q < 3, 



(2 - 9)P y^^ftF (fEf 



\y\ 



2y/(q-l)ft 

\y\ 



and for — 1 < q < 1, 



A, 



2(2 - g)/? ^2(1 



2v/(l-?)/9 



(20) 



(21) 



with i/ = ^-4 — h- J-v is the Bessel functions of first kind and K u is the 
modified Bessel function of second kind. 

It is possible to have variations of the g-mexican hat by using ft = /3(g) 
(with ft(l) = 1), for instance, ft = 1/(3 — g). 



3 Modulated (/-Gaussian 

Now we turn to the Morlet's wavelet, or modulated Gaussian [JUl EDI OH] > a 
function associated with the birth of wavelet analysis 15^] . Within the 
context we are dealing with, we search for a wavelet that is a generalization 
of HOI EH: 

h{x) = 7T" 1 / 4 (e~ lkx - e" fc2 / 2 ) e~ x2 '\ (22) 

where k = vr(2/ ln.2) 1 / 2 . So we simply modulate the usual trigonometric 
functions (e~ thqX ) with a g-Gaussian (e"^ ): 

h q (x) = B q (e~ ik " x - A q (k q )) e"^ 2 , 00 < q < 3. (23) 

The function A q (k q ) is such that the admissibility condition, written in the 
domain of frequencies, 

F[h q (x); 0] = 0, (24) 
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is satisfied. It means that 



AJk q ) = 



T 


e 


^g 


; o 


T 


y* 2 ; o 





(25) 



Taking into account the Fourier transform of a g-Gaussian (see Equations 
3.384 9 and 3.387 2 of HSj), and T[(T ik « x ] y] = V2tt 5(y - k q ), and also, 
from the convolution theorem (with the symmetric convention we are 
adopting here), 



F[f(x)g(x);y] = F(y)*G{y) 

= -L F(y - G(0 d£, 

V27T J-oo 



we find for q > 1, 



(26) 



and for g < 1, 



A g (fc 9 



1 - g 2 



2^(1-9)/? 



fen 



J-, 



(27) 



(28) 



with /i = — |. We follow the same criterion adopted by 40 a for the 
determination of the value of k q . the ratio between the second highest and 
highest local maxima of Re h q is |. It results that 



(29) 



The normalization constant B q is given by 



.9-1 2/ 



1/2 



(i- g ) 1/2 r(iVf) 



1/2 



1 < q < 3, 
g < 1. 



(30) 



(B q is an approximation, as B\ = 7r -1 / 4 of Equation (|22|) is also an approxi- 
mation.) Figure 121 illustrates v^Re h q with (3(q) = 1/(3 — q). (The factor \[2 
is used to have the real part normalized.) We observe the long (power-law) 
tail for 1 < q < 3, in marked contrast with the rapidly vanishing tail for 
q — 1. The cut-off is present for q < 1. In the case q — > — oo, h q (x) reduces 
to a two-cycle function, the imaginary part of it is a variation of the one-cycle 
sine presented in |4*D*] . 



4 g- Trigonometric wavelets 

The g-exponential function Q, expanded to the imaginary (or, more gen- 
erally, complex) domain by analytic continuation, leads to g-trigonometric 
functions 32 : 



COSg X 



e g e q 



2i 



The g-cosine and g-sine functions may be expressed as e 
where 



p q {x)e 



At?) 



3 (i-?)^ 2 



arctani[(l — q)x] 



(31) 

itp q (x) 

(32) 



We want to construct a wavelet based on such g-trigonometric functions. 
For this purpose, we recall that the derivative of a g-exponential may be 
expressed as 



de| 
dx 



e 2-l/<7 • 



(33) 



Once p q> i(x) — > for \x\ — > oo, the admissibility condition is satisfied for 1 < 
g < 2. Renaming the parameter g, we define the following g-trigonometric 
wavelet: 



wt q (x) = C q Bq 



1 < q < 2. 



(34) 



The normalization constant is given by: 



a, 



i i 



2-q Tr 1 ^ 



(*-i)r A 



r ( — — - 

1 q-l 2 



n V2 



(35) 
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For brevity of notation, we can write the real and imaginary parts of wt q (x) 
as 



wc q {x / 



WS q {X 



= V2 Re wt q (x) = \p2C q cos g ' ( 36 ) 

1 = V2 Im wt q {x) = V2C g sin, (kZT^) ' ( 37 ) 



Some typical curves are shown in Figure E2 (for brevity, we only show the 
odd function ws g (x)). Note that the number of oscillations decreases as q goes 
from 1 to 2. The functions present infinite oscillations of vanishing amplitudes 
at q — > 1 (C g _»i — > 0). wSg(x) presents only one root, at xo = 0, for | < q < 2, 
and wCq(x) presents only one pair of roots, at x = ± tan [(g — , 
for | < g < 2. As g — > 2, C g — > 00, and the wavelets become sharply 
localized (the roots of wc q (x) approaches ±2/ it). Let us also mention that 
the modulation of the functions is not Gaussian, but rather a power-law 
decay. But they are essentially different from the modulated g-Gaussian 
derived in the previous Section. 

The Fourier transform of wc q (x) and ws q (x) are found with Equation (j31j) 
and (see Equations 3.382 6 and 3.382 7 of 

y>o, 

2vr {-y)^~ l n (38) 

y < 0, 




[(?-i)a]^r(^) 



and 



2n yi- 1 e^ -1 ' 11 



-, y>o, 



^foT ;y] = ^ [(g-iJal^r^)' ' (39) 

0, y<0. 



5 An Example 

Wavelet transforms constitute a multi-purpose tool that rapidly met wide- 
spread applications in many areas of pure and applied sciences. To exemplary 
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demonstrate the reliability of g-wavelets, we briefly present how they can be 
used to reproduce several properties of well known fractal sets. 

Consider a self-affme profile y = f*(x) shown in Figure 0] Its construction 
starts with the generator shown in the Inset. The generator is recursively ap- 
plied for each straight line segment. The self-affme fractal profile is obtained 
in the limit of infinite iterations. 

The resulting figure after infinite generations has a local fractal dimension 
D = 2 — a, where a = \nB y / \nB x = 0.5 is the roughness exponent, and 
B x and B y are the scaling factors along the x and y axis respectively. To 
analyze the scaling properties of f(x) around an arbitrary point xo, we shift 
the origin and define 

U {x) = f{x + x)- /(x ) « Ax<*» (40) 

where a(xo), the Holder (or singularity) exponent of f Xo {x), indicates how 
this function vanishes (or diverges) at x = xq- 

Using a general wavelet transform of f Xo (x), it is straightforward to show 
that the Holder exponent follows immediately from the wavelet transform as 

T g (a,b) f X0 (x)^a a ^\ (41) 

Figure El shows how the q- wavelet transform (Equation (}4~Tj) ) behaves for 
two points along the profile for some values of q. The curves clearly indicates 
that the value a = 0.5 is accurately reproduced. We observe that the results 
for q > 1 show small irregularities. This is due to the fact that, as e~ x decays 
only as a power-law when x —>■ oo, the numerical integration in (|41|) must 
be performed over a wider interval than for q < 1. In order to make evident 
this effect, we have performed all integrals over a fixed integration interval. 
On the other hand, results for q < 1 show good convergence, as the function 
naturally presents a cutoff. Oscillations in the curves are typical for this kind 
of analysis, as shown in the same figure for the usual mexican hat (q = 1). 

Results for a more complex situation can be explored if we consider a 
multifractal set. It is generated along a similar way used to obtain the first 
set, where we choose two different scaling factors for the first and second half 
of the profile in each generation of its construction. The plots for different 
values of xq have different slopes (Figure EJ), indicating that many scaling 
laws and fractal dimensions are found in the resulting profile. 

q- Wavelets also prove to be a reliable tool to analyze this much more 
complex situation that requires a multi-fractal formalism. This amounts to 
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evaluate the generalized fractal dimensions Dq, or its Legendre transformed 
singularity spectrum f{a) = r(q) — aQ, where a represents the Holder 
exponent. This analysis proceeds within the so-called Wavelet Transform 
Modulus-Maxima Method (WTMMM), that has been developed in recent 
years. We will not show the details that can be found, e.g., in |55j . 

Figure [7| shows the f(ce) spectra, for different values of q, corresponding 
to the multifractal profile. The graphs indicate that the spectra produced 
by q- wavelets, with q < 1, are comparable to the one obtained by the usual 
q = 1 mexican hat. As we fixed the integration interval for all values of q, we 
observe major differences in the spectra for values q > 1. This limitation to its 
practical use is expressed by larger numerical effort to compute the integrals 
in the wavelet transforms with the same accuracy for values of q < 1. 

6 Final Remarks 

The functions here introduced generalize widely used mother wavelet func- 
tions by means of a single parameter. The g-wavelets present significant 
different behaviours, when compared with the original q = 1 cases. The 
fingerprints of the generalized q- wavelets are the power-law decay (q > 1), 
and the cut-off (q < 1). The generalization was inspired in the nonextensive 
statistical mechanics and the g-exponential that emerges from it. Within 
this formalism, the entropic index q measures departures from the usual 
Boltzmann-Gibbs behaviour. There is a number of possible origins for such 
departures, like slow (power-law) mixing, long-range (spatial) interactions, 
long-term (temporal) memory or fractal nature of the phase space. When 
dealing with signals which exhibit some of these features, one can hopefully 
take advantage of the index q to tune the wavelet to the signal, as it is here 
briefly exemplified for the q- mexican hat. Other wavelets, as well as scaling 
functions, may be generalized along the lines here introduced, for instance, 
the Meyer wavelet or the sine function. 
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Figure 1: g-mexican hat with (3 = 1/2. (a) Normalization constant A q ; (b) 
ip q (x) for 1 < q < 3; (c) < q < 1; (d) — 1 < g < 0. The usual mexican hat 
(g = 1) is represented by a dotted line, for comparison. 
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Figure 2: Normalized real part of the g-modulated Gaussian (hq(x) = 
V^Re h q ). (a) q = 1 (usual case); (b) q = 2; (c) q — > — oo (illustrated 
with q = —100). Note that abscissa scale in (c) is different from the others. 
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Figure 5: Wavelet transform (Equation (|4ip) with two different values of b: 
x l Q = 0.25 and x r Q = 0.50. Note that the curves for the same values of q 
coincide, as it should be for a monofractal. 
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Figure 6: Wavelet transform (Equation (|41|) with two different values of b: 
x l = 0.25 (upper curves) and x r = 0.75 (lower curves). The slopes are clearly 
different, as a signature of the multifractal nature. 
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Figure 7: f{ct) spectra for the multifractal profile. For q = 0.8 results are 
very similar to those with q — 1. q = 0.6 yields spurious results. 
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